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Abstract. We describe the implementation of weak quantum measurements in su¬ 
perconducting qubits, focusing specifically on transmon type devices in the circuit 
quantum electrodynamics architecture. To access this regime, the readout cavity 
is probed with on average a single microwave photon. Such low-level signals are 
detected using near quantum-noise-limited superconducting parametric amplifiers. 
Weak measurements yield partial information about the quantum state, and corre¬ 
spondingly do not completely project the qubit into an eigenstate. As such, we use 
the measurement record to either sequentially reconstruct the quantum state at a 
given time, yielding a quantum trajectory, or to close a direct quantum feedback 
loop, stabilizing Rabi oscillations indefinitely. 


Measurement-based feedback routines are commonplace in modern elec¬ 
tronics, including the thermostats regulating the temperature in our homes to 
sophisticated motion stabilization hardware needed for the autonomous oper¬ 
ation of aircraft. The basic elements present in such classical feedback loops 
include a sensor element which provides information to a controller, which in 
turn steers the system of interest toward a desired target. In this paradigm, 
the act of sensing itself does not a priori perturb the system in a significant 
way. Moreover, extracting more or less information during the measurement 
process does not factor into the control algorithm. For quantum coherent cir¬ 
cuits, these basic assumptions do not hold. The act of measurement is invasive, 
and the so-called backaction drives a system toward an eigenstate of the mea¬ 
surement operator. Furthermore, the information content of the measurement 
determines the degree of backaction imparted. For example, strong measure¬ 
ments completely project a superposition state into a define eigenstate while 
extracting enough information to unambiguously allow an observer to deter¬ 
mine which eigenstate has been populated. Weak measurements on the other 
hand, accrue less information, indicating which eigenstate is more likely to 
be populated, and have proportionally weaker backaction that does not com- 
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pletely collapse a quantum superposition. This type of measurement provides 
a natural route to implement active feedback in quantum systems with the 
weak measurement outcome providing the input for a quantum controller that 
takes into account the measurement induced backaction. 

In this chapter, we first discuss a general formalism based on positive 
operator-valued measures (POVMs) which can describe both weak and strong 
measurements. We then apply this formalism to a superconducting two-level 
system coupled to a microwave frequency cavity, which is well approximated 
by the Jaynes-Cummings model commonly used in cavity quantum electro¬ 
dynamics. In such a system, the measurement strength for a given integration 
time can be adjusted by varying the number of photons in the cavity. To access 
the weak measurement regime, the cavity is typically populated, on average, 
with less than one photon, requiring an ultra-low-noise amplifier for efficient 
detection. In contemporary experiments, this function is realized using su¬ 
perconducting parametric amplifiers. After briefly describing these devices, 
we discuss two basic types of experiments. In the first set, the result of a 
sequence of weak measurements is used to reconstruct individual quantum 
trajectories using a Bayesian update procedure. The statistical distribution of 
an ensemble of many such trajectories is then analyzed to experimentally and 
theoretically obtain the most likely path. In these experiments, the backaction 
results either solely from the measurement process or from the combination 
of measurement and unitary evolution under a coherent drive. In principle, 
the reconstructed state and the detector values corresponding to the most 
likely path can be used for arbitrary feedback protocols and optimal control. 
In the second type of experiment, we use the weak measurement outcome to 
directly complete an analog feedback circuit. In particular, we demonstrate 
the real-time stabilization of Rabi oscillations resulting in the suppression of 
their ensemble decay. Finally, we close with a discussion of future directions 
for hardware improvement and additional experimentation, particularly in 
multi-qubit systems. 


1.1 Generalized Measurements 

In quantum mechanics, predictions about the outcome of experiments are 
given by Born’s rule which for a state vector | ipj) provides the probability 
P(a ) = \(a\ipi}\ 2 that a measurement of an observable described by an oper¬ 
ator A with eigenstates |a) yields one of the eigenvalues a. As a consequence 
of the measurement, the quantum state is projected into the state |a). Here 
we consider a qubit, with states |0) and 11), which is conveniently described 
by the Pauli matrix algebra with pseudo-spin operators a x , a y , and a z . For 
example, if we prepare the qubit in an initial state, 

i^) = i+^)^4(i°)+i 1 )) 


(i.i) 
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Born’s rule tells us the probability of finding the qubit in the state |0), 
P(+z) = 1/2. Here we have introduced the notation ±x, ±z to indicate the 
eigenstate of the a x and a z pseudo-spin operators. After the measurement, 
the qubit will remain in the eigenstate corresponding to the eigenvalue and 
subsequent measurements will yield the same result. In this way the mea¬ 
surement changes ( a z ) from 0 to 1, while changing ( o x ) from 1 to 0. This 
change, associated with a projective measurement of a z is the backaction of 
measurement and also referred to as the collapse of the wavefunction. 

The fact that the measurement takes (a x ) from 1 to 0 is a consequence of 
the Heisenberg uncertainty principle. Because a x and a z do not commute, if 
one component of the pseudo-spin is known then the others must be maxi¬ 
mally uncertain. Thus some amount of backaction occurs in any measurement, 
and measurements (such as the projective measurement discussed above) that 
cause no more backaction than the amount mandated by the Heisenberg un¬ 
certainty principle are said to be quantum non-demolition. 

In this chapter, we are concerned with a more general type of measurement 
that can be performed on the system [1-6]. Real measurements take place over 
a finite amount of time and we are interested in describing the evolution of 
the system along the way. As such, a projective measurement is composed 
of a sequence of partial measurements. Formally, we can describe these mea¬ 
surements by the theory of positive operator-valued measures [1-4], (POVMs) 
which yields the probability P(m ) = Tr(12 m pf2' m ) for outcome m, and the as¬ 
sociated backaction on the quantum state, p P(m) for a system 

described by a density matrix p. POVMs are “positive” simply because they 
describe outcomes with positive probabilities and “operator-valued” because 
they are expressed as operators. These operators fl m obey = I 

as is necessary for POVMs. When fi a = |a)(a| is a projection operator and 
p = \ip)(il)\, the theory of POVMs reproduces Born’s rule. For the case de¬ 
scribed above, the projective measurement of a z is described by the POVM 
operators Q +z = ( I+a z )/2 and 1?_ 2 = (I—a z )/2 and straightforward applica¬ 
tion of the theory reproduces the effects of projective measurement discussed 
above. 


1.1.1 Indirect Measurements 

While it is tempting to connect a quantum system directly to a classical mea¬ 
suring apparatus, such an arrangement often results in far more backaction 
on the quantum system than is dictated by the Heisenberg uncertainty princi¬ 
ple. For example, an avalanche photodiode can be used to detect the presence 
of a single photon, but it achieves this detection by absorbing the photon, 
thereby completely destroying the quantum state. In practice, classical de¬ 
vices are composed of too many noisy degrees of freedom to perform QND 
measurements. 

To circumvent this problem, we break the measurement apparatus into 
three parts to conduct an indirect measurement. The three parts of the appa- 
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ratus are the quantum system of interest, a quantum pointer system that can 
be coupled to the quantum system, and then a classical measurement appara¬ 
tus that can be used to record the pointer system. To execute a measurement, 
the pointer system is coupled to the quantum system and after sufficient in¬ 
teraction the two systems become entangled. Then, a classical measurement 
apparatus, such as a photodetector or an ammeter, is used to make measure¬ 
ments on the pointer system. These measurements on the pointer system are 
projective and depending on their outcome, the quantum system’s state is 
accordingly changed. 

To illustrate this process consider a system consisting of two qubits [3]. We 
initialize one qubit (the system qubit) in the state in (1.1) and allow a second 
(environment) qubit to interact with the primary qubit for a brief period of 
time such that the two qubits are in the entangled state, 

W) OC ((1 + e )|0)sys + (1 — e)|l)sys) ® |0) e nv 

+ ((1 ~ c)|0)sys + (1 + e)|l)sys) ® |l)env- (1-2) 

We then make a projective measurement of the second qubit (using the POVM 
l?±z), and if the result is the |0) state, then the system qubit is left in the 
state, 


W) a (1 + e)|0) sys + (1 — e)|l) sys . (1.3) 

For e <C 1 this partial measurement slightly drives the system qubit to the 
ground state. However, what happens if we instead measure the environment 
qubit in the a y basis, |y±) = ^(|0) ± z11)? We can express the joint state in 
this basis, 

|<F) oc ((1 + e - i(l - e))|0) sys + (1 - e - *(1 + e))|l) sys ) (g> |y+) e nv 

+ ((1 — £ — i{ 1 + e))|0) sys + (1 + € — *( 1 — e))|l)sys) ® |j/-)env- (1-4) 

By factoring out a global phase factors <f> + = 7r/4 + e and </>_ = 7t/4 — e and 
assuming that eCl, the joint state is, 

|<F) oc e* 0+ (|0) sys + e _2 * e |l) sys ) ® |y + ) e nv 

+ (|0)sys + e 2?£ |l) sys ) ® |y-)env (1-5) 

which shows that the system qubit state obtains a slight rotation depending 
on the detected state of the environment qubit. This simple model indicates 
two salient features of partial measurements. First, if the environment and 
the system are weakly entangled, then projective measurements on the envi¬ 
ronment cause only small changes in the qubit state, and second, the choice 
of measurement basis for the environment gives different conditional evolu¬ 
tion for the system. After the projective measurement on the environment 
qubit the system is in a product state with the environment qubit and no 
entanglement remains. 
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We now imagine many such successive interactions between the system 
qubit and a set of environment qubits. In this case each entangling interaction 
and projective measurement of the environment causes a small amount of 
random backaction on the qubit, which results in a diffusive trajectory for the 
state of the qubit. Such weak measurements have been recently implemented 
using coupled superconducting qubits [7]. 

1.1.2 Continuous Measurement 

The measurements that are realized in the cQED architecture described later 
in section 1.2 can be represented by a set of POVMs {fly}, where V is the 
dimensionless measurement result which is scaled so that it takes on average 
values ±1 for the qubit in states ±z respectively, 

f2 v = (2na 2 )- 114 e^ v ~^ 2 / 4a2 \ ( 1 . 6 ) 

Here, l/4a 2 = hr]At, and k parametrizes the strength of the measurement, p 
is the quantum efficiency, and At is the duration of the measurement. The 
operators fly satisfy f fly fly dV = I as expected for POVMs. The probability 
of each measurement yielding a value V is P(V) = Tr (fly p t fly), which is 
the sum of two Gaussian distributions with variance a 2 centered at +1 and 
— 1 and weighted by the populations poo and pn of the two qubit states. 
The a z term in fly causes the back action on the qubit degree of freedom, 
p —> flypfly/P(V), due to the readout of the measurement result V. Eq. 
(1.6) can also describe a stronger measurement, ultimately yielding the limit 
where the two Gaussian distributions are disjoint, and the readout projects 
the qubit onto one of its a z eigenstates, with probabilities poo and pn. The 
strength of the measurement is controlled by the parameter k and we will see 
later in section 1.2 how this quantity is related to experimental parameters. 

For weak measurements, At <C t c , where the characteristic measurement 
time t c = l/4/c?y describes the time it takes to separate the Gaussian measure¬ 
ment histograms by two standard deviations. The distribution of measurement 
results is then approximately given by a single Gaussian that is centered on 
the expectation value of a z , 

P(fly) cz e~ 4kriAt{v - {rT ^ )2 , (1.7) 

which highlights the fact that V is simply a noisy estimate of ( a z ). This 
allows the measurement signal to be written as the sum of ( a z ) and a zero 
mean Gaussian random variable [2]. 

The time evolution of the quantum state following a sequence of measure¬ 
ments described by Eq. (1.6) in the limit At —> 0 is given by the stochastic 
master equation [1]: 

^ = k(a z pa z - p) + 2rjk(a z p + pa z - 2Tr (a z p)p)V ( t ). (1.8) 
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Fig. 1.1. Circuit QED setup consisting of a qubit dispersivley coupled to a readout 
cavity. A signal transmitting the cavity at the cavity frequency aqcuires a qubit- 
state-dependent phase shift, shown as two phasors in the I — Q plane. Phase sensitive 
amplification amplifies one quadrature exclusivley. 


Here, the first term is the standard master equation in Linclblad form and 
the second term is the stochastic term that updates the state based on the 
measurement result. In the case of unit detector efficiency (?y = 1) the stochas¬ 
tic master equation (perhaps surprisingly) yields pure state dynamics. This 
happens because the decoherence term with rate k in the equation is exactly 
compensated by the stochastic term to yield a random pure state evolution; 
the first term shrinks the Bloch vector towards the z-axis, while the stochastic 
term adds a random transverse component, putting it back to the surface of 
the Bloch sphere again. In contrast, if we set rj = 0, all random backaction is 
suppressed, and we instead obtain a conventional, deterministic master equa¬ 
tion, where probing of the system causes extra Lindblad decoherence on the 
system, but yields no information. 


1.2 Quantum Measurements in the cQED Architecture 

The circuit QED architecture [8] provides both excellent coherence properties 
and high fidelity measurement. In particular, it is an ideal test bed for study¬ 
ing quantum measurements because it enables one to implement text book 
quantum measurements relatively easily. The two key aspects which make this 
possible is the applicability of an ideal measurement Hamiltonian [8] and the 
availability of quantum limited parametric amplifiers [9,10]. 

Circuit QED (cQED) is essentially the implementation of the cavity QED 
architecture [11] using superconducting circuits. Instead of an atom inter¬ 
acting with the electromagnetic field inside a Fabry-Perot cavity, the cQED 
architecture uses a superconducting qubit as an ‘artificial’ atom which inter¬ 
acts coherently with the electromagnetic field in an on-chip waveguide res¬ 
onator [12] or a 3D waveguide cavity [13]. Figure 1.1 shows the basic cQED 
setup. Dispersive readout in cQED was introduced in Chapter 6 and we briefly 
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reintroduce it here to make this chapter more self contained. The Hamiltonian 
describing this interaction is the Jaynes-Cummings Hamiltonian, 

H = —^ -<j z + TiLo c (a) a + -) + hg(aa + + aV_) (1.9) 

where the first term represents the qubit as a pseudo-spin, the second term 
is the cavity mode, and the third term represents the interaction between the 
qubit and the electromagnetic field in the rotating wave approximation. Here 
woi is the transition frequency between the qubit levels, w c is the cavity mode 
frequency, g is the coupling strength between the qubit and cavity mode and 
<j± = ( a x ±ia y )/2 are the qubit raising and lowering operators. Typically, the 
qubit frequency is far detuned from the cavity frequency to protect the qubit 
from decaying into the mode to which the cavity is strongly coupled. In this 
dispersive regime (A = woi — lo c g), the effective Hamiltonian reduces to, 

H = + h{u c + xcr z )(a)a + (1.10) 

where \ is called the dispersive shift. It is now possible to see how the mea¬ 
surement is implemented. The cavity mode frequency is a function of the 
qubit state with y, which depends on g and A, setting the magnitude of this 
shift. The measurement proceeds by probing the cavity with a microwave 
signal and detecting the qubit state dependent phase shift of the scattered 
microwave signal. In this architecture, one can clearly see the process of in¬ 
direct measurement introduced earlier. The pointer system is the microwave 
field which interacts with the qubit and gets entangled with it. One then mea¬ 
sures the output microwave field which in turn determines the qubit state. As 
mentioned earlier, the exact measurement on the pointer determines the back- 
action on the qubit and that can be controlled by the choice of amplification 
method used to detect the output microwave field. We will however come to 
that a little later. 

1.2.1 Dispersive Measurements 

We will now formalize the interaction between the qubit and the pointer state 
using the language of coherent states to describe the quantum state of the 
microwave field. Let |a) represent the coherent state of the microwave field 
sent to the cavity with an average photon number n = |a-| 2 . The initial state 
of the qubit is oo|0) + ai|l). After interaction, the final entangled state of the 
system is given by 

IV’/) = d 0 \e~ l9d a)\0) + ai\e +l9d a)\l) (1.11) 

where 9d is the dispersive phase shift of the scattered microwave signal. The 
two coherent states can be represented in the quadrature plane as shown in 
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Fig. 1.2. Output coherent states represented in the IQ plane. The length of the 
vector increases with increasing photon number in the coherent state and leads to 
better distinguishability between the qubit states. 

figure 1.2. For a fixed dispersive phase shift 6d, the larger the average number 
of photons in the coherent states, the easier it is to distinguish them. Here, 
the distinguishability of the output coherent states is directly related to the 
distinguishability of the qubit state i.e. the strength of the measurement can 
be controlled by the microwave power. 

One can now choose to measure the output microwave field by performing 
a homodyne measurement and choosing one of the two possible quadratures. 
From figure 1.2, it is clear that the qubit state information is encoded in the 
‘Q’ quadrature while the T quadrature has the same value for both qubit 
states. We will first describe the single quadrature measurement of the ‘Q’ 
quadrature using a phase-sensitive amplifier and the associated backaction on 
the qubit. Since the coherent state is not a two-level state, a measurement on 
it yields a continuous range of values and corresponds to the case discussed 
in section 1.1.2. Typically, experiments with microwave signals employ mixers 
operating at room temperature to implement homodyne detection. However, 
the microwave signals used to probe the cavity are extremely weak and need 
to be amplified before they can be processed by room temperature electronics. 
The challenge is that commercial amplifiers add significant noise which results 
in imperfect correlation between the measured output and the measurement 
backaction on the qubit. We use superconducting phase-sensitive parametric 
amplifiers to implement near noiseless amplification of a single quadrature of 
the microwave signal. 

1.2.2 Parametric Amplification 

There are several parametric amplifier designs [9,10,14] but we will restrict 
our discussion to one based on a lumped element non-linear oscillator [10]. The 
basic physics of parametric amplification in such a system can be understood 
by considering a driven, damped Duffing oscillator model [15] whose classical 
equation of motion is given by 



( 1 . 12 ) 
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Fig. 1.3. (a) The basic circuit diagram of a Josephson junction based non-linear 
oscillator. A circulator is used to separate the amplified reflected signal from the 
incident signal, (b) Phase diagram of a driven non-linear oscillator. The solid line 
indicates the effective resonant frequency as a function of drive power. The dashed 
lines enclose the bistable regime of the non-linear oscillator. The oval marks the 
paramp biasing regime, (c) The reflected pump signal phase is plotted as a function 
of pump power and defines the transfer function of the paramp. 


Here ujq is the linear resonant frequency for small oscillations, r is the ampli¬ 
tude damping coefficient, 5 is the gauge-invariant phase difference across the 
junction which is the dynamical variable of oscillator, and cod is the frequency 
of the harmonic driving term often referred to as the pump. 

Figure 1.3(a) shows the basic circuit diagram of a Josephson junction 
based non-linear oscillator. The phase diagram of such a non-linear oscillator 
as a function of drive frequency and drive power is shown in Figure 1.3(b). 
A characteristic feature of a non-linear oscillator is that its effective resonant 
frequency is no longer independent of the amplitude of its oscillations. This 
is shown schematically by the black solid line in figure 1.3(b) where the effec¬ 
tive resonant frequency decreases with increasing driving power for Josephson 
junction based nonlinear oscillators. Beyond a critical drive frequency and 
power (w c and P c ), the driven response becomes bistable and we are not in¬ 
terested in that regime. The relevant regime for parametric amplification is 
marked in Figure 1.3(b) and is just before the onset of bistability. Because 
of the power dependence of the effective resonant frequency, it is possible 
to cross the resonance (black line) both in frequency and amplitude. Figure 
1.3(c) shows the phase of the reflected pump signal from the oscillator as a 
function of pump amplitude. This is essentially the transfer function of the 
amplifier and the pump power is chosen to bias the system in the steep part 
of this curve. 

The signal to be amplified has the same frequency ujd with an amplitude 
which is typically less than 1% of the pump amplitude, is combined with 
the pump signal and sent to the nonlinear resonator. Since the signal and 
pump are at the same frequency, one can define a phase difference between 
them. If the signal is in phase with the pump, then the combined amplitude 
changes significantly with the input signal which moves the bias point and 
consequently the reflected pump phase. This is essentially the mechanism of 
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amplification, i.e. a small signal brings about a large phase shift in the large 
pump signal. However, if the signal is 90 degrees out of phase with the pump, 
then to first order, it has no effect on the bias point and consequently no 
change in the reflected pump phase. So the signals which are in phase with 
the pump get amplified while the quadrature phase signals get de-amplified. 
This allows one to selectively amplify the ‘Q’ quadrature of the microwaves 
scattered from the cavity in a cQED measurement. In principle, such a phase- 
sensitive amplifier can amplify one quadrature without adding any additional 
noise to the signal [ 10 ], resulting in a signal-to-noise-ratio that is maintained 
after amplification and an output noise level that is entirely set by amplified 
zero point fluctuations associated with the coherent state. In practice, due to 
signal losses between the cavity and the amplifier and due to inefficiencies in 
the amplifier itself, one typically obtains an efficiency r] ~ 0.5. In other words, 
the output noise is only about double the unavoidable quantum noise [16,17]. 

1.2.3 Weak Measurement and Backaction 

A single weak measurement V m is obtained by integrating the homodyne signal 
Vq (t) corresponding to the ‘Q’ quadrature for a measurement time r, 

Vm ^ = ^yJ 0 V Q^ dt (!- 13 ) 

where AV is the difference in the mean value of Vq (t) corresponding to the 
two qubit states. This implies that the mean value of V m is ±1 for the two 
qubit states. Measuring the ‘Q’ quadrature of the scattered microwave signal 
in cQED using a phase-sensitive amplifier corresponds to a cr z measurement of 
the qubit. The corresponding backaction pushes the qubit state towards one 
of the eigenstates of a z . This is shown as the solid black line in figure 1.4(a) 
where Z z = (& z )\v m is plotted as a function of weak measurement result V m 
for an initial qubit state (|0) + \\))/\/2. Similary, the expectation values of 
a x and a y are also shown as solid blue and red lines respectively. The dashed 
line corresponds to the theoretical prediction given by [17,18], 

= tanh j = \J 1 - (Z z ) 2 er' yT , Y z = 0 (1.14) 

where r is the measurement duration, S = 64rx 2 n?7/K is the dimensionless 
measurement strength and 7 is the environmental dephasing rate of the qubit. 
To connect with the nomenclature used in section 1.1.2, we note that V m 
is equal to the dimensional measurement result V defined earlier whereas 
S = 4 /a 2 where a 2 is the variance of Gaussian measurement distributions. If 
77 = 1, the state remains pure during the entire measurement process. Here, 
the data corresponds to rj = 0.49 and S = 3.15. 

It is important to note that even though the other quadrature ‘I’ does 
not contain any information about a z , it does contain information about the 



1 Quantum Trajectories 


11 



Fig. 1.4. Backaction of a weak measurement, (a) When the ‘Q’ quadrature of 
the scattered microwave signal is amplified a a z measurement of the qubit occurs. 
The dashed lines indicate the theoretical prediction from (1.14) and the solid lines 
are the tomographic results conditioned on the measurement result V m . (b) A ‘<f>’ 
measurement of the qubit is obtained by amplifying the T quadrature. The dashed 
lines indicate the predictions from (1.15) and the solid lines are the conditioned 
tomography. The data correspond to r/ = 0.49 and S = 3.15 and correspond to the 
experimental setup in [17]. 


photon number fluctuations in the coherent state [18]. This implies that a 
measurement of the ‘I’ quadrature will result in the qubit state rotating about 
the Z axis in the Bloch sphere [17] and will not push the state towards one of 
the eigenstates of a z . This can be seen in figure 1.4(b) where the expectation 
values of cr x ,y, z are plotted (solid lines) after a weak measurement with the 
paramp amplifying the T quadrature. Here V m is obtained by integrating the 
homodyne signal Vj(t) corresponding to the ‘I’ quadrature for a measurement 
time r. Note that this is not the same as measuring the a x or the a y operator 
and we label it as a “</>’ measurement since it results in rotation about the Z 
axis in Bloch sphere. Consequently, there is no asymptotic value for or 
Y* and they evolve sinusoidally as, 

X* = cos ^pp^j e _7r , Y+ = - sin e" 7T , Z* = 0. (1.15) 

and Z$ doesn’t evolve and remains zero. 

A measurement using a phase-sensitive amplifier truly enables one to mea¬ 
sure one quadrature only while erasing the information in the other quadrature 
so that no observer can get access to that information. In fact, that informa¬ 
tion no longer exists and is not the same as ignoring the information in that 
quadrature which would lead to decoherence. So a single quadrature measure¬ 
ment of the ‘Q’ quadrature only provides information about a z and there are 
no photon number fluctuations. Similarly, when measuring the ‘I’ quadrature, 
no information about cr z is available and hence (<r z ) does not change [17]. 

It is also possible to measure both quadratures simultaneously by using 
a phase preserving amplifier. The parametric amplifier described above can 
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Fig. 1.5. A sequence of weak measurements leading to a quantum trajectory of 
the qubit state on the Bloch sphere. The qubit is initially along x = 1 and each 
measurement imparts conditional dynamics on the state. 




be used in phase preserving mode by detuning the signal frequency from the 
pump frequency [10] or one could use a two mode amplifier like the Josephson 
Parametric Converter (JPC) [14]. In this case, each measurement gives two 
outputs corresponding to the T and ‘Q’ quadrature [19]. This implies that one 
learns about the qubit state as well as the photon number fluctuations and 
consequently the measurement backaction results in the Bloch vector rotating 
around the Z axis while it approaches one of the eigenstates of a z . Somewhat 
surprisingly, even in this case the qubit state can remain pure [19] throughout 
the measurement process provided the measurement efficiency r) = 1. 


1.3 Quantum Trajectories 

Quantum trajectories were first introduced as a theoretical tool to study open 
quantum systems [20-24]. Rather than describe an open quantum system by 
a density matrix, which for a TV-dimensional Hilbert space requires TV 2 real 
numbers and requires solving the master equation, the evolution of a pure 
state (which requires only TV complex numbers), can be repeatedly calculated 
to determine the evolution of p(t). The quantum trajectory formalism thus 
assumes that the evolution of an open (and therefore mixed) quantum sys¬ 
tem can be expressed as the evolution of several, individual, pure quantum 
trajectories. 

1.3.1 Continuous Quantum Measurement 

The process of continuous quantum measurement can be built out of a se¬ 
quence of discrete weak measurements as sketched in figure 1.5. As we have 
already established, each of these measurements induces a specific conditional 
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Fig. 1.6. Individual quantum trajectories of the qubit state as given by z = {a z ) 
and x = (a x ) from an initial state +x. Four different trajectories are shown in color 
on top of a greyscale histogram indicates the relative occurrence of different states 
at different times. 


backaction on the quantum state. If several of these weak measurements are 
conducted in series then the discrete time evolution of the quantum state can 
be determined. The cQED measurement apparatus that we consider forms a 
continuous probe of the quantum system, however since the cavity has a finite 
bandwidth re, the measurement signal is correlated at times less than ~ 1/re 
and it makes sense to bin the measurement signal in time steps At ~ 1/re. 
For a cavity with k/2tt = 10 MHz, this correlation time is roughly 16 ns, so 
the continuous measurement record is discretized in similar time steps. These 
discrete, weak measurements can be considered to be continuous in the limit 
where the discretization steps are much smaller than the characteristic mea¬ 
surement time t c = K/16nx 2 nr) which was defined in section 1.1.2. Typically, 
experiments operate with r c ~ 1 /rs, such that At <C r c . 

A single experimental iteration results in a continuous measurement sig¬ 
nal V (t) which is subsequently binned into a string of measurement results 
(Vi, Vi+i, I 1 + 2 ). Given the initial state of the qubit, the state is updated at 
each time step U, either based on the stochastic master equation [2,25], or a 
Bayesian argument [17]. This leads to a discrete time trajectory x(t),y(t), z(t). 
Figure 1.6 displays several of these trajectories for the qubit initialized in the 
state x = +1. Because each experiment results in a different measurement 
signal V(t), each trajectory is different. 

To verify that these trajectories, which are conditional on a single mea¬ 
surement signal, are correct we have to prove that at every point along the 
trajectory the conditional state makes correct predictions for the outcome of 
any measurements that can be performed on the system. To accomplish this, 
we perform conditional quantum state tomography at discrete times along the 
trajectory. We denote a single trajectory x(t), z(t) as a “target” trajectory. 
This trajectory makes predictions for the mean values of measurements of a x 
and cr z versus time. To create a tomographic validation of these predictions at 
a specific time U we perform several experiments with a weak measurement 
duration of ti that conclude with one of three tomographic measurements. 
For each of these experiments, we calculate the trajectory x'(ti) , z' (ti) and 
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Fig. 1.7. Tomographic reconstruction of the trajectory. To reconstruct the target 
trajectory at time z(ti) several experiments are performed with a projective mea¬ 
surement of a z at time ti. For each of these, the state (x'(ti ), z'(t)) is calculated 
and if it matches the target trajectory the outcome of the projective measurement 
is included in the average. 


if x'{ti) = x(U ) ± e and z'[ti) = z(U) ± e the outcomes of the projective 
tomographic measurements are included in the average trajectory x(t),z(t). 
We find that the target trajectory and tomographic trajectory are in close 
agreement. 

If we did not condition our trajectories on the measurement signal, for 
example if ij = 0, then the trajectories would simply follow the ensemble 
evolution as described by a standard Lindblad master equation. In this case, 
each trajectory would be the same, yet would still make correct predictions 
for the outcome of projective measurements performed on the system. While 
this unconditioned evolution makes correct predictions, it also quickly takes an 
initial pure state to a mixed state. However, if ?y ~ 1, the conditional quantum 
state retains substantial purity for all time and makes correct predictions for 
measurements on the system. 

1.3.2 Unitary Evolution 

So far in this section we have considered the case of quantum non-demolition 
(QND) measurement, in that the weak measurements we perform oc <j z com¬ 
mute with the qubit evolution Hamiltonian. Since the measurements are QND, 
all the measurements commute and only the integrated measurement signal is 
necessary for the conditional state evolution. In this regime, the measurements 
can be treated simply in terms of classical probabilities, since the evolution 
only depends on the state populations, and the measurement signal oc <j 2 also 
depends only on the populations. 

The situation becomes much more rich if we allow for unitary evolution 
of the qubit state that does not commute with the measurement operators. 
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Fig. 1.8. Quantum trajectories in the presence of unitary driven evolution. Indi¬ 
vidual quantum trajectories of the qubit state as given by z = ( <j z ) from an initial 
state +x. Six different trajectories are shown in color on top of a greyscale histogram 
indicates the relative occurrence of different states at different times. 


To accomplish this, we drive the qubit resonantly to induce Rabi oscillations. 
This drive is described by the Hamiltonian, Hu = hf]<j y /2. In this case, the 
unitary evolution is included in the stochastic master equation, 

^ = -^r[#R,p] + k{a z pa z - p) + 2pk(a z p + pa z - 2Tr (a z p)p)V t . (1.16) 

The Rabi drive turns coherences into populations and vice versa, causing the 
measurement signal to depend on the coherences of the qubit. Such evolution 
is fully quantum and reveals interesting features regarding the competition 
between unitary dynamics and measurement dynamics. Figure 1.8 displays 
several trajectories which exhibit this competition between measurement and 
driven evolution. The trajectories are oscillatory, but distorted by the stochas¬ 
tic backaction of the measurement, approaching jump-like behavior expected 
for the regime of quantum jumps. This figure also highlights how state tracking 
maintains the purity of the state in comparison to the full ensemble evolution. 


1.3.3 The Statistics of Quantum Trajectories 

We have so far demonstrated the ability to track individual quantum tra¬ 
jectories which evolve in response to measurement dynamics (wave function 
collapse) and also in competition with unitary driven evolution. These tra¬ 
jectories are stochastic in the sense that the evolution of each trajectory is 
different. But, what statements can we make about these trajectories in gen¬ 
eral? Clearly, the evolution associated with weak measurement would be in 
some way different than the dynamics of quantum jumps [16], but how do we 
quantify and characterize trajectories? 

In order to characterize the general properties of trajectories we need a 
method to look at ensembles of trajectories but that does so in a way that 
depends on the individual trajectories. To accomplish this, we consider the 
sub-ensemble of trajectories that start and end at certain points in quan¬ 
tum phase space. This sub-ensemble is pre- and post-selected in that all the 
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trajectories start from the same initial state (pre-selection) and we then post¬ 
select a sub-ensemble that ends in a specific final state. Given these pre- and 
post-selected trajectories we can now examine aspects of the sub-ensemble. 
One such property is the most-likely path that connects the initial and fi¬ 
nal states. This particular choice is of interest because such most-likely paths 
can be calculated with a stochastic action principle for continuous quantum 
measurement which maximizes the total path probability connecting quantum 
states [26]. Experiments [27] show good agreement with the theoretical most- 
likely path and the predicted path from theory, thus validating the theory 
which may be applicable in other quantum control problems. This analysis 
gives insight into the dynamics associated with open quantum systems, with 
applications in quantum control and state and parameter estimation. 


1.3.4 Time-symmetric State Estimation 

The examination of pre- and post-selected quantum trajectories raises the 
notion of time symmetry in quantum evolution and quantum measurement 
[28-31]. We have so far used the quantum state as predictive tool, that is, at a 
time t the quantum state described by p(t) makes correct predictions for the 
probabilities of the outcomes of measurements performed at time t and the 
associated mean values of these observables. However, after a measurement is 
performed, the quantum state may continue to evolve due to further probing 
and unitary driving, and we may ask at some later time T > t what is the 
probability for the outcome of that measurement in the past given the results 
of later probing. 

Consider a measurement scenario where two experimenters monitor the 
evolution of a qubit and track its quantum state p(t). At time t one experi¬ 
menter makes a measurement of the qubit but locks the result “in a safe”. The 
second experimenter then continues to monitor the qubit and at a later time 
T the second experimenter wants to guess the outcome of the measurement 
whose result is locked in the safe. Clearly more information is available if the 
second experimenter accounts for information about the qubit obtained after 
the first measurement, and if this experimenter can correctly account for those 
results, he will be able to make more confident predictions for the result in the 
safe. Stated simply, the second experimenter must determine what result is 
most likely to be locked in the safe given the subsequent measurement signal. 

One can show [32] that at time T, the second experimenter’s hindsight 
prediction for the measurement performed in the past is given by, 

P p (m) = Tr (MmP(t)f2lE(t)) 

£ m Tr (n m p(t)rf m E(t)) 

which describes the probability of obtaining outcome m from the POVM 
measurement fl m performed by the first experimenter. Here p{t) is the usual 
quantum state propagated forward in time until time t and E(t) is a similar 
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matrix which is propagated backwards from time T to time t using a sim¬ 
ilar method for the calculation of p(t). The matrix E(t) has recently been 
calculated for experimental data and demonstrated that Eq.(1.17) makes cor¬ 
rect and indeed more confident predictions for the outcome of measurements 
performed in the past [25,33]. 


1.4 Analog Feedback Stabilization: Rabi Oscillations 

In this section, we will explore how one can use the continuous measurement 
record obtained using weak measurements to modify the behavior of the quan¬ 
tum system being monitored with the help of feedback. Unlike feedback in 
classical systems, one has to worry about the random backaction associated 
with the measurement of a quantum system. However, as explained in previ¬ 
ous sections, the state of a quantum system can be monitored perfectly using 
weak measurements if the initial state is known and the measurement effi¬ 
ciency is unity. Even though the evolution of the quantum system is random 
and unpredictable, it is still knowable by monitoring the weak measurement 
output. This can also be done in the presence of any additional unitary evo¬ 
lution as explained in 1.3.2. However, decoherence processes which are always 
present will tend to take the system away from a desired state one might 
want. We will now describe how one can use quantum feedback to prevent 
the quantum system from deviating from a desired state. We will consider the 
particular case of feedback control in a resonantly driven qubit undergoing 
Rabi oscillations [34]. 

1.4.1 Weak Monitoring of Rabi Oscillations 

Let us first look at Rabi oscillations more carefully. The state of a resonantly 
driven qubit evolves sinusoidally between its two states with a rate 12# which 
depends on the strength of the resonant drive. To be specific, for a qubit state 
a(t)|0) +/3(t)|l), |a(f)| 2 = sin 2 (12#f/2). These oscillations in the qubit state 
probability are called Rabi oscillations. In the absence of any decoherence, 
given the initial state and the Rabi frequency 12#, we can predict the qubit 
state at any future time. One can equivalently say that the phase of the Rabi 
oscillations is known and remains unchanged with time. However, decoherence 
processes will introduce errors in this deterministic evolution and over some 
characteristic time scale, the phase of the Rabi oscillations will diffuse. 

In a typical Rabi oscillation experiment, the qubit is initialized in the 
ground state and resonantly driven for a fixed duration of time (r#) followed 
by a projective measurement. This process is repeated many times to obtain 
the ensemble averaged qubit state (a z ). A plot of (cr z ) as a function of r# 
yields decaying sinusoidal oscillations where the decay constant depends on 
the decoherence in the system. The decaying oscillations are an indication 
of the diffusion of the Rabi oscillation phase with time. Here, qubit driving 



18 


K. W. Murch, R. Vijay, and I. Siddiqi 


and projective measurement are never done simultaneously. However, one can 
drive the qubit while measuring it weakly and still obtain ensemble averaged 
Rabi oscillations. As discussed in 1.1.2, since the weak measurement output 
can be thought of as a noisy estimate of (<r z ), an ensemble average of the weak 
measurement signal in the presence of Rabi drive also yields decaying oscilla¬ 
tions. Figure 1.9 shows Rabi oscillations obtained using weak measurements. 
However, one important difference is that the decay constant now depends on 
both environmental decoherence and measurement strength. This additional 
measurement induced decoherence is a consequence of the ensemble averag¬ 
ing where we ignore the individual results of the weak measurements. This is 
in contrast to the oscillatory quantum trajectories shown in Figure 1.8 with 
simultaneous Rabi driving and weak measurement. 

We will now discuss a feedback protocol [34] which corrects for the phase 
diffusion of Rabi oscillations and prevents the decay of Rabi oscillations. In 
principle, one can do a full reconstruction of the quantum state [35] to estimate 
the feedback signals required. However, to do that in real-time is experimen¬ 
tally challenging. Instead, we use classical intuition in this feedback protocol 
motivated by phase-locked loops (PLL) used to stabilize classical oscillators. 
In a PLL, one compares the phase of an oscillator with that of a reference 
signal. Any phase error is then corrected by creating a feedback signal propor¬ 
tional to the error which controls the oscillator frequency. Essentially, if the 
oscillator is lagging in phase, then the feedback signal increases the frequency 
and vice-versa. 

We can now apply the same idea to our weakly monitored qubit in the 
presence of Rabi driving [16]. The basic feedback setup is shown in Figure 
1.10. A reference signal at the Rabi frequency is multiplied with the weak 
measurement signal and low pass filtered to create the error signal. Since the 
weak measurement signal is a noisy oscillatory signal corresponding to the 
oscillating qubit state, the error signal is proportional to the deviation in 
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Fig. 1.9. Rabi oscillations obtained using ensemble averaged weak measurements. 
The decay constant is set by both environmental decoherence and measurement 
induced decoherence. 
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Fig. 1.10. Simplified feedback setup for stabilization of Rabi oscillations. 
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Fig. 1.11. Feedback stabilized Rabi oscillations. Initially the ensemble averaged 
Rabi oscillations decay, primarily due to measurement induced decoherence. Feed¬ 
back is turned on later and results in the oscillations recovering as the phase of Rabi 
oscillations synchronize with the reference signal. 



phase of the Rabi oscillations with respect to the reference signal. The error 
signal is used to control the amplitude of the Rabi drive which in turn controls 
the Rabi frequency just as in a PLL. Figure 1.11 shows the effect of such a 
feedback signal which is turned on after a time much greater than the decay 
constant of the Rabi oscillation in the absence of feedback. On can clearly see 
that the ensemble averaged oscillations recover when the feedback is turned on 
and stabilize to a fixed amplitude. As long as the feedback is on, the ensemble 
averaged oscillations will never decay. This implies that the Rabi oscillations 
have synchronized with the reference signal and the phase diffusion has been 
reduced due to feedback though not completely eliminated. Note that the 
slow drift in the mean level of the signal is due to finite probability of getting 
excited into the second excited state of the transmon qubit [16]. 

To ensure that the stabilized oscillations are not an artifact of the mea¬ 
surement setup, quantum state tomography was used to verify the quantum 
nature of the stabilized oscillations. Figure 1.12(a) shows a plot of (ax), (ay), 
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Fig. 1.12. (a) Quantum state tomography of feedback stabilized Rabi oscillations 
showing (ux), (&y), and (rrz) for one full Rabi oscillation, (b) Feedback efficiency 
is plotted a function of feedback strength. 


and ( az) for one full stabilized Rabi oscillation. The magnitude of these oscil¬ 
lations do not reach ±1 indicating that the synchronization is not perfect and 
the phase diffusion of the Rabi oscillations has not been completely eliminated. 
The data shown is the best synchronization we obtained in this experiment 
corresponding to a feedback efficiency D = 0.45 which is approximately given 
by the amplitude of the oscillations in (crz) or (crx)- The best synchroniza¬ 
tion is obtained for an optimal feedback strength F as is evident from Figure 
1.12(b) which shows a plot of D vs F (solid squares). The dimensionless feed¬ 
back strength F is essentially the feedback loop gain and is controlled by the 
amplitude of the reference signal. 

There are two main factors in this experiment which result in D < 1. The 
first one is the total measurement efficiency ijt which had a value of 0.4 in this 
experiment. This efficiency has two contributions given by r/ t = r] rj env . Here 77 
is the measurement efficiency due to the detector as introduced in 1.1.2 and is 
set by the excess noise introduced by the amplification chain. The other term 
Venv = (1 + T , env /T , m ) _1 takes into account the environmental dephasing T env . 
where r m is the measurement induced dephasing. In order for feedback stabi¬ 
lization to work well, one needs to ensure that the environmental dephasing is 
small compared to the measurement induced dephasing i.e. r] e nv — > 1. In other 
words, you want the measurement induced disturbance to dominate over the 
environmental disturbance since the measurement output can then be used 
to correct for those. The second factor affecting the feedback efficiency is the 
loop delay. Since the qubits are operating inside a dilution refrigerator while 
the feedback electronics are operating at room temperature, there is a delay 
in the feedback signal which results in inefficient synchronization. The solid 
line in Figure 1.12(b) is obtained using numerical simulations including the ef¬ 
fect of feedback delay and shows good agreement with the experimental data. 
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While one might be tempted to increase measurement strength arbitrarily to 
approach ry t = 1, feedback delay and finite feedback bandwidth leads to an op¬ 
timal value of measurement strength for maximizing feedback efficiency. This 
is because stronger measurement requires faster feedback which is limited by 
the bandwidth of the feedback loop. 

This experiment demonstrates the use of continuous measurement and 
feedback to stabilize Rabi oscillations in a qubit. The simple feedback proto¬ 
col which is based on classical intuition works successfully because the feed¬ 
back signal achieves near perfect cancellation of the random measurement 
backaction for an optimal value of F. This technology can provide another 
route for quantum error correction based on weak continuous measurements 
and feedback [36, 37] with a potential advantage in situations where strong 
measurements can cause qubit state mixing [38]. Such techniques also offer 
the possibility of measurement based quantum control for solid-state quantum 
information processing [39-45]. 


1.5 Conclusion 

Weak measurements realize a flexible method of implementing active feedback 
in quantum systems. With superconducting circuits that operate at microwave 
frequencies, feedback fidelity is currently limited by the overall measurement 
efficiency of the amplification chain. In typical setups, the amplifier is housed 
in a separate cryo-package and inefficiencies result from losses in cable con¬ 
nectors and passive components such as circulators which add directionality 
to the signal path. One promising avenue to overcome this limitation is to 
use parametric devices that can be directly integrated on-chip with the qubit. 
Such types of circuits use a combination of complex pumping techniques at 
different frequencies [46] or multiple cavity modes to isolate the amplifier bias 
form the qubit circuit [47]. Additionally, as more sophisticated feedback rou¬ 
tines are developed, particularly sequences which involve digital processing of 
the measurement data, then loop delays resulting from long data paths and 
latencies in classical electronics must also be taken into account. Wiring com¬ 
plexity and dead time maybe minimized by integrating quantum circuits with 
cryogenic classical logic, either superconducting or semiconducting. 

Another frontier to be explored in superconducting circuits is the use of 
feedback and control in multi-qubit arrays. For example, in such an archi¬ 
tecture, one can imagine simultaneously performing weak measurements on 
each qubit. The resulting joint-state information can potentially be used to 
reconstruct the initial state of the array from families of quantum trajecto¬ 
ries, realizing another form of state tomography that capitalizes on the dense 
information embedded in the correlations of an analog weak measurement. 
Such techniques can also be extended for Hamiltonian parameter estimation, 
parity measurements, and measurement based error correction. In essence, si- 
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multaneous probing of a quantum many body system parallelizes both ‘read’ 
and ‘write’ operations. 
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